Magnetic breakdown induced Peierls transition 
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We predict the new type of phase transition in quasi one-dimensional system of interacting elec- 
trons at high magnetic fields, the stabilization of a density wave which transforms a two dimensional 
open Fermi surface into a periodic chain of large pockets with small distances between them. We 
show that quantum tunneling of electrons between the neighboring closed orbits enveloping these 
pockets transforms the electron spectrum into a set of extremely narrow energy bands and gaps 
that decreases the total electron energy, thus leading to a magnetic breakdown induced density wave 
ground state analogous to the well-known instability of Peierls type. 
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Interacting electrons with quasi one-dimensional 
(Q1D), i. e. open and corrugated, Fermi surfaces show in 
many cases instabilities leading to a density wave (DW) 
ground state with a broken translational symmetry [1]. 
This type of phase transition, predicted by Peierls in 1955 
[2], occurs due to divergent charge fluctuations which 
open the gap at the Fermi surface and so lower the en- 
ergy of the new ground state. While the dominance of 
a repulsive Coulomb or an attractive phonon mediated 
electron-electron interaction are as a rule responsible for 
respective spin or charge modulation of DW, its stabi- 
lization depends on the degree of nesting between the 
left and right two dimensional Fermi surfaces shown in 
Fig. 1. Nesting in real materials is always imperfect, 
leading to two-dimensional [3] pockets of finite size after 
determining the optimal wave vector of DW ordering, 
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Q = (2kp + 5k\\, b*/2 + 5kx), where k F and b* are Fermi 
wave number and transverse reciprocal lattice vector re- 
spectively, while <%i_i_ are small corrections. Unlike the 
nested parts of the Fermi surface, these pockets develop 
a weaker gap in the electron spectrum, or even remain 
gapless. The pockets thus act against the DW, and, if 
large enough, may completely eliminate it. 

A qualitative change however takes place after apply- 
ing strong enough magnetic field H perpendicular to the 
plane (p x ,p y ) in Fig. 1. Landau quantization of band 
states then reduces the transverse derealization of elec- 
trons, making the system "more one-dimensional" [4] and 
therefore more susceptible to the DW ordering. DWs are 
then strengthened, or even newly established in systems 
which otherwise remain metallic down to T — 0. The lat- 
ter case of field-induced density wave (FIDW) is partic- 
ularly relevant for our further considerations. The best 
known examples are those realized in a series of Bech- 
gaard salts with spin modulation [5] , in which the FIDW 
order takes place in available magnetic fields up to ~ 50T, 
provided the deviation form perfect nesting is not large. 
The pockets with discrete Landau orbits are then small 
and, after establishing DW gap in dominant "nesting" 
parts of Fermi surface, quite distant one from another. 



FIG. 1: Quasi one-dimensional Fermi surface, with arrows 
denoting directions of semiclassical electron motion under a 
magnetic field perpendicular to the plane (p x ,p y ). b% and by 
are reciprocal lattice constants. 



The Landau quantization leads to the localization of 
electron motion, therefore increasing the band energy in 
general, or leaving it unchanged for the discrete set of 
values of H for which the Fermi energy lays exactly in 
the middle between two neighboring Landau levels [6]. 
Still, it contributes to the FIDW stabilization through an 
additional gain in the correlation energy part due to the 
improved longitudinal localization of electron states. On 
the other hand, the magnetic breakdown effect [7, 8], i. 
e. the magnetic field assisted electron tunneling through 
barriers between neighboring pockets which dclocalizes 
the electron motion and so decreases the band energy, is 
far too weak to be relevant. 

At this point we come to our central question, namely, 
is it possible to stabilize a density wave order with the 
help of orbital quantization in the quasi one-dimensional 
band, but through a predominant decrease of the band 
energy due to the magnetic breakdown? This magnetic 
breakdown induced density wave (MBIDW) is possible 
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FIG. 2: (a) Periodic net of large closed orbits separated by 
small classically inaccessible MB arrears (thick dots). Q — 
(Q ~ 2A;_f,0) is lattice deformation wave vector, and the dis- 
tance between closed orbits is proportional to the product of 
electron mass and deformation potential V(x). (b) Areas of 
closed orbits, 5i,2, and areas under open orbits, S3, 4 for the 
respective left and right motion of electrons, determine the 
phases of semi-classical electron wave functions defined in the 
text. 

only by strengthening the tunneling between neighboring 
pockets, i. e. by choosing the DW wave vector which, 
oppositely to FIDW in Bechgaard salts, favors short and 
weak tunnel barriers together with large pockets between 
them. In the present work we examine this new possi- 
bility by using the standard mean field approach for DW 
instability [9] and the semiclasical treatment of MB tun- 
neling [8], and establish the conditions under which this 
gain overwhelms the clastic energy loss due to the lattice 
deformation [10], leading so to the MBIDW stabilization. 

Qualitative considerations. Let us choose the defor- 
mation wave vector Q = (Q w 2kp,0) creating a static 
periodic lattice distortion and a charge modulation V(x) 
with the period 2-kTi/Q that combines open trajectories 
with opposite directions of the electron motion, shown by 
arrows in Fig. 2. We encounter a chain of closed orbits 
encircling large pockets and small areas between them 
(shown with thick dots) at which MB takes place. Here 
we treat MB electron dynamics within the semi-classical 
approximation: an electron moving along a semi-classical 
section of the chain, say 1, is scattered at the MB area to 
section 2 and section 2' with the amplitude probabilities p 
and t respectively, fulfilling the normalization condition 
Ip| 2 + l r | 2 = 1 W1 th the MB probability given by 

H 2 =cxp(- 7 )- (1) 

IV I 2 

where 7 = 7r— — Here v x and v v are longitudinal 
and perpendicular projections of the electron velocity 
v = de(p)/dp at the point of the magnetic breakdown, 
a = TieH/c, e and c are the electron charge and the light 
velocity, respectively, while Vq is the non-diagonal matrix 



FIG. 3: Energy spectrum for electrons shown in Fig. 2 in the 
regime of strong MB (|r| 2 < 1). n is the MB band number 
and P x o is the conserved generalized electron momentum. The 
period of E n (P x o) is 2%a / B y v, while the widths of the energy 
bands and gaps are AE^and ~ p| 2 ^c and AE gap ~ |r| 2 fio; c 
respectively (uj c is the cyclotron frequency). Dotted straight 
lines denote the electron spectrum in the absence of MB. 

element of V(x). For the electron-phonon system within 
the mean field approximation, \Vo\ = 2gb where g is the 
coupling constant and b =< 5q > is the mean value of 
the phonon annihilation operator at the wave number Q 
(see e.g. [9]), here directly proportional to the amplitude 
of the lattice displacement. Thus, due to the finite lat- 
tice distortion V(x) electron scatters in a periodic set of 
regions of MB, latter having a role analogous to that of 
atomic lattice in an one-dimensional metal. 

In order to understand qualitatively the structure of 
the energy spectrum of an electron on the orbit chain let 
us assume at first that |V(:r)| is negligibly small while 
the strength of the magnetic field is finite. In this case 
the probability of MB is \p\ 2 = 1 and the electron moves 
along an open trajectory 1-2-3 ... or l'-2'-3'... (trajec- 
tories I and II, respectively). It is known [11] that the 
electron spectrum E^' II (P x0 ) of such motion is continu- 
ous and characterized by a discrete quantum number n 
and a continuous momentum P x q. Here superscripts refer 
to trajectories I and II; P x q is the generalized momentum 
projection conserved in the gauge with the vector poten- 
tial A = (— Hy, 0,0). The average quantum-mechanical 
electron velocities are v 1 ' 11 = dE^' 11 (P x0 ) / dP x0 (v 1 > 
and v 11 < 0). This electron spectrum is represented by 
the dotted straight lines in Fig. 3. 

For a finite charge modulation V(x), the MB is not 
total (\p\ 2 < 1) and the trajectories I and II are coupled 
by the finite tunneling probability. The degeneracy at 
the crossing points in the electron spectrum is then lifted 
by |t| 2 , producing there energy gaps analogous to those 
in the spectrum of a real one-dimensional electron gas 
created by the charge modulation V(x) itself [12]. 

If the Fermi energy occurs in the middle of one of these 
gaps, that is at one of the crossing points in the un- 
perturbed spectrum, the electron band energy decreases. 
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This energy gain grows by increasing the gap, which in 
turn increases with |r|. As |r| increases by increasing \ Vq\ 
as is seen from Eq. 1, the energy band gain is favored by 
the increase of the charge modulation V(x). However, by 
increasing V(x) one increases the lattice energy and, as a 
result, competition of these two tendencies, a decrease of 
the band energy and an increase of the elastic one, gives 
the optimal value of the charge modulation V(x). 

Analytical calculations. In the p-space the electron 
moves along classical trajectories p x \p y ) (marked by 
I = 1, 1 , ... in Fig. 2) between the points of MB. Corre- 
sponding semi-classical wave functions are 
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where p x \p y ) are solutions of the equation s s {p x ,p y ) ~ 
E while e s {p) is the electron dispersion law for a finite 
V(x) and at H = 0, s = 1,2 indexing the new elec- 
tron bands. At the points of magnetic breakdown the 
in-coming and out-coming functions (2) are coupled by 
the magnetic breakdown 2x2 unitary matrix [8, 13, 14] 



-IT) 



(3) 



e. g. a 2 = exp(i^)(ai \p\ exp(w) + ay |r| exp(ir])). The 
phases depend on the MB parameter 7 from Eq. 1. In 
the limit 7^1 (that is |t| 2 <C 1) to which we limit our 
calculations of MBIDW, 77 f« -7r/4-7ln7, v « 0,* « 0. 

Matching wave functions (2) at the MB points with 
matrix (3) and taking into account the electron phase 
gains between them one finds the dispersion function 

P no* 

D(E,P x0 ) = cos$_ - |r| 2 cos$+ - |p| 2 cos($ - -^). 

a 

(4) 

Here 2*, $± = (l/2a){S 1 (E) ± S 2 (E)) and 

$0 = (l/2a)(S 3 (E)-S 4 (E))+2v, while S h2>3A (E,Q) are 
the areas enclosed by electron trajectories which depend 
on the position of the crossing points of trajectories I and 
II determined by Q (see Fig. 2). The electron spectrum 
is determined by dispersion equation D(E, P x q) = 0. 

If \p\ 2 = 0, this equation gives a set of discrete Landau 
levels, while in the case \p\ 2 7^ the degeneracy with 
respect to P x o is lifted and the electron spectrum depends 
on a continuous quantum number, E — E„(P x q) 7 where 
the band number n enumerates solutions of the dispersion 
equation. Solving this equation in the limit |r| 2 <C 1 one 
finds the electron spectrum presented in Fig. 3. 

In order to find an explicit formula for the density of 
states (DOS) for electrons moving along the MB chain, 
we express it in terms of the spectral function (4) as fol- 
lows: 

f dP 
vmb{E)= jJ25(E-E n (P x0 ))^> 



\D (E,P x0 )\6(D(E,P xQ )) 
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dP x0 
2itTi ' 
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where here and below / = df/dE. Integrating with 
respect to P x q we find 
a 
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where 0(x) equals unity for x > and zero otherwise. 
From here one easily sees that at |r| 2 = (that is at Vq = 
0) DOS is equal to the one in the absence of MBIDW . 

Eq. 6 permits to find the number of electrons N = 
J[exp(E — fi)/T + \]~ 1 v{E)dE and the electron free en- 
ergy F = Np, -T f ln[l + exp(^ - E) /T]v(E)dE, n being 
the chemical potential. We calculate them expanding the 
right-hand side of Eq. 6 in a double Fourier series in <f>-|- 
and in a power series in |r| 2 . Below we solve the problem 
in the vicinity of the critical temperature T c in which the 
order parameter and the potential Vb are small, so that 
|t| <C 1 and we can keep only terms of the lowest order 
(~ |t| 2 ) in the above-mentioned power series. 

After finding TV and F we take into account the conser- 
vation of electron number and determine the change of 
the chemical potential Sp = (Imb — MOj where \imb and 
Ho are chemical potentials at V(x) ^ and V(x) = 
respectively. Inserting it in the free energy and neglect- 
ing terms of the order of \t\ a Tilu c / £f and T 2 (huj c /eF) 2 , 
we find the difference of the total free energies in the 
presence and the absence of MBIDW at T > Tilj^ /it, 
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where wJ T ' = eH/m^c with the effective mass given by 
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Here S± = Si(p 07 Q) ± 5 2 (/Jcb Q) is the effective electron 
loop area, and the last term is the lattice elastic energy 
given by phonon frequency uiq at momentum Q. Eq. 7 
is valid in the regime a < \S-\ <C 5+. 

From Eq. 7 it follows that AF has a series of minima 
at the discrete set of lattice deformation wave numbers 
Q n for which S-(no, Q)/2a = (7r/2)(2n+ 1) with integer 
values of 77 (for dSi/dQ w -dS 2 /dQ) [15]. On the other 
hand, the minimization of Eq. 7 with respect to b 2 shows 
that A_F is negative and has an absolute minimum at 
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as shown in Fig. 4. Here A = g 2 / (huiQ n po) is the dimcn- 
sionlcss electron-phonon coupling constant (see [9]). The 
corresponding critical temperature for the transition to 
the MBIDW state is 



fi^j lnf l C0S (' S '+/ 2cr ) 
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FIG. 4: The free energy difference AF/AFq vs. x 5 
(ng 2 /a\v x v y \)b 2 for 7 = (a\v x v y \/ng 2 )(hujQ / AF ) = 0.8 
showing stable minimum of DW ground state at xq w 0.22. 



provided A/zo |cos(S , + /2cr)| /m^lu^l > 1. 

In conclusion, wc have shown that the MBIDW insta- 
bility may occur in quasi-one-dimcnsional conductors in 
a wide range of parameters and at any finite magnetic 
field. The topological prerequisite for such instability 
is the open Fermi surface of the type shown in Fig. 1, 
which after a periodic breaking of symmetry from Fig. 
2, leads to the lattice of closed orbits connected by bar- 
riers at cross points. It was already pointed out that 
magnetic field assisted tunneling through these barriers 
is crucial for the stabilization of MBIDW. MB introduces 
a qualitative change in electron dynamics, transforming 
the continuous spectrum of quasi one-dimensional inter- 
acting electrons under a strong magnetic field into a set 
of alternating narrow energy bands and energy gaps, with 
the widths proportional to cyclotron energy multiplied by 
electron MB tunneling and scattering probability respec- 
tively. Just a mere Landau quantization of such closed 
electron orbits does not lead to the DW stabilization since 
it only rises the band energy due to the localization of 
electron motion. Only by derealization of motion intro- 
duced by MB, whose negligence was probably the reason 
why this effect was not predicted in preceding literature, 
the total energy is lowered enough to grant the stability 
of DW ground state. 

We have also shown that the free energy of MBIDWs 
has local minima at a set of discrete values of the defor- 
mation wave vector Q, indicating the role of resonances 
between neighboring orbits with generally different size 
(Fig. 2). We remind that this series of wave vectors 
belong to the "anti-nesting" regime, which again reflects 
the essential differences between MBIDWs and standard 
FIDWs, the latter being the consequence of logarithmic 
anomalies in the DW correlation function due to the al- 



most perfect nesting between left and right Fermi surfaces 
in Fig. 1. The dependence of T c on H has a periodic 
(in 1/H) set of gaps for a given Q n (Eq. 9), in con- 
trast to the corresponding dependence for FIDWs [16]. 
This, together with a qualitatively different dependence 
on the coupling constant and band parameters, could be 
one guide in the search for MBIDWs. Another one are 
fast magneto-resistance oscillations due the formation of 
relatively large closed orbits (Fig. 2). The regime wider 
than that of Eq. 7 and the peculiarities of MBIDW phase 
diagram will be analyzed in a more detailed presentation. 
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